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PREFACE 



This final report is the summary of a two-year effort portraying 
the development, calibration and verification of a one-dimensional variable 
cross- sectional area numerical model. This is one of the publications 
of the thermal pollution group at the University of Miami. 

This effort covers the applications of this model to Cayuga Lake 
in New York and Lake Keowee in South Carolina. The conclusions and 
derivations are accompanied by an abundance of figures and tables. 
These figures and tables are presented in such a way as to make it 
possible for use In calibrating or verifying other one-dimensional models. 

This work was conducted under funding from National Aeronautics 
and Space Administration (NASA) , Kenedy Space Center and the Environ- 
mental Protection Agency (EPA). 



ABSTRACT 



A one-dimensional model for studying the thermal dynamics of cool- 
ing lakes has been developed and verified. The model is essentially a 
set of partial differential equations which are solved by finite difference 
methods. The model includes the effects of variation of cross- sectional 
area with depth, surface heating due to solar radiation absorbed at the 
upper layer and internal heating due to the transmission of solar radia- 
tion to the sub-surface layers. The exchange of mechanical energy be- 
tween the lake and the atmosphere is included through the coupling of 
thermal diffusivity and wind speed. The effects of discharge and intake 
by power plants are also included. 

The numerical model was calibrated by applying it to Cayuga Lake. 
The performance was very good. The model; was then verified through 
a long term simulation using Lake Keowee data base. The comparison 
between measured and predicted vertical temperature profiles for the 
nine years is good. The physical limnology of Lake Keowee is presented 
through a set of graphical representations of the measured data base. 
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SECTION 1 
INTRODUCTION 



Deep bodies of water provide a convenient source of condenser 
cooling water supply to electric generating power plants. One of the 
problem, associated with this, that has gained increasing importance in 
recent vears is the thermal pollution caused by the discharge of waste heat 
from these power plants. The degradation of the quality of these waters 
usually occurs through the direct influence of the increased temperature 
on aquatic life or through the lowering of the amount of dissolved oxygen 
and other biochemical effects. 

In temperate regions most deep bodies of water develop a thermo- 
cline during their annual heating cycle. A warmer epilimnion at the top 
is isolated from a cooler hypolimnion below by severe stable gradients 
(thermocline) . Convective transport and heat addition caused by power 
plant discharge result in disturbances in the thermocline. 

The formation time, phasing, depth and severity of the thermocline 
are crucial factors affecting the bio-chemical processes in an aquatic 
ecosystem. The nutrient levels, species spectra and physical character- 
istics are quite different in the two distinct domains below and above the 
thermocline. This stratification is believed to be caused by nonlinear 
Interaction between the wind- generated turbulence and stable buoyancy 
gradients. While being heated from above, a basin forms stable strati- 
fication thereby inhibiting wind-generated turbulence. During early 
scring, most temperate lakes exhibit a nearly homothermal temperature 
distribution with a temperature of about 4°C (which is the temperature 
of maximum density for water) extending all the way to the bottom. As 
the air above the lake begins to warm, the lake receives heat, at an 
increasingly rapid rate. During the early part of the warming season, 
the lake continues to remain nearly homothermal, since the heat that is 
received at the surface layers is transported to the deeper layers by 
wind-induced currents and turbulence. As the rate of heating of the 
lake continues to increase, the rate at which heat is received at the 
surface layers soon exceeds the rate of heat removal to the deeper lay- 
ers, and the temperature of the surface layers begins to increase. 
During this period the temperature decreases monotonically with increas- 
ing depth. As heating continues, a point of inflexion develops in the 
temperature profile separating an upper well-mixed layer from a lower 
well-mixed layer. The region around the point of inflexion is a region 
of intense stable temperature gradient and consequently, low turbulence 
levels. Turbulent diffusion through this region is minimal. Further 



heating merely increases tha temperature of the well-mixed upper layer. 

As cooling begins, the wind mixing is augmented by convective 
mixing caused by static instability and the thermocline (the region con- 
taining the point of inflexion) recedes into the deeper layers of the lake. 
The epilimnion cools until the lake reaches its minimum heat content and 
near homot hernial conditions result. 

This report describes the application of a one-dimensional model in 
pr<idicting the above phenomena. This model was calibrated using data 
from Cayuga Lake and the calibrated model was applied to Lake Keowee. 
The purpose of the effort is to develop a predictive ability for long-term 
thermal behavior of cooling lakes. 



SECTION 2 
CONCLUSIONS 



A one-dimensiona; model which includes area change with depth, 
vertical convection, varying diffusivity, thermal discharges, and internal 
absorption of radiation has been developed. The model was calibrated 
using Cayuga Lake data base and verified with Lake Keowee data base. 

The application to Cayuga Lake indicates excellent performance. 
A comparison of cylindrical and paraboloid approximations for the lake 
indicate significant differences in thermocline depth, eddy-diffusivity, 
and temperatures at mid-depths. This shows that effects of area change 
with depth are not negligible. The insignificant surface temperature 
difference is attributable to the fact that equal surface areas were used 
in both simulations. However, these effects will be more pronounced 
in real basins where decrease in area with depth is more severe than 
the linear variation for the paraboloid case. 

The long term (nine-year) simulations of temperature profiles, 
formation and decay of the thermocline for Lake Keowee compare fairly 
well with measured data. For these simulations, the model did not need 
recalibration. While the influence of discharge and ambient conditions 
were satisfactorily modeled, tie calculated and measured profiles had 
much larger differences for Lake Keowee than for the one-year simulation 
in Cayuga Lake. Two reasons are believed to cause this relatively great- 
er inaccuracy. The first is the complicated discharge patterns in Lake 
Keowee caused by power plant, hydro and pumped storage operation 
compared to ambient simulation for Cayuga Lake. The second cause is 
that for Cayuga Lake; the assumption of horizontal mixing is more mean- 
ingful since it is a single basin domain, while Lake Keowee is a double 
basin domain connected by a canal. The model, however, can be adapted 
for application to two connected basins. 



SECTION 3 
RECOMMENDATIONS 



The main disadvantage of a one- dimensional thermal model ties in 
the fact that resolution is sacrificed for computational speed. Three di- 
mensional models are bulky and time consuming but have much better 
resolution, however, when long term simulations are necessary, a one- 
dimensional model is recommended. 

The model described here can be modified to include the single 
effects of the various quantities involved in the surface heat transfer 
phenomenon rather than using the equilibrium temperature concept. 
This is particularly recommended for the user who is interested in 
modeling the long term effects of one (for example, evaporation) of 
the quantities involved in the surface heat transfer processes. 

Furthermore, the model can be easily adapted to handle connected 
multiple domains. This recommendation is discussed in the text. 



SECTION 4 
THE MATHEMATICAL MODEL 



DESCRIPTION 

The one-dimensional variable cross- sectional area model used for 
this study is described in this section. The detailed derivation of the 
governing equations is described in Appendix A, while the boundary 
conditions and the important terms occurring in the governing equations 
are discussed here. 

This model which assumes lateral uniformity was developed to pre- 
dict the vertical temperature profiles for complete annual cycles in a lake. 
The model includes the effects of variation of the horizontal cross-section- 
al area with depth. Surface heating due to solar radiation absorbed at 
the surface layer and the internal heating due to the transmission of the 
unabsorbed solar radiation to the deeper layers of the lake are taken 
into account. The exchange of mechanical energy between the- lake and 
the atmosphere is accounted for through the friction velocity and the 
thermal diffusivity, since these quantities can be related to wind speed. 
Finally, the effects of power plant discharge and intake are considered. 

The full three-dimensional equations of mass and heat balance are: 
and 



n = -7-o7 (U.I) 



h ^ C « T) s 7-oC R-7T - 7-pC TV + H (1.2) 

3t p p p 

Applying the divergence theorem, and then differentiating and integrat- 
ing, a set of one-dimensional equations are obtained: 

AU) || = | 1 Mz)pV z + IA' (4.3) 

and 

A(z) |- t (pCpT) =|- z fpC p A(r)K z |X] - |~ z (pC A(i)TV z ) +QA' + A U)H (z) 

(4.4) 
where, 

z is the vertical coordinate, measured upward from the deepest point of 



of the lake. As a subscript it marks the vertical component of a 
vector . 

A(z) is the horizontal cross-section of the lake at height z. 

p is the density. 

t is the time. 

V is the vertical velocity of flow. 

C is the heat capacity. 

T is the temperature . 

K is the vertical component of 'K\ the heat diffusivity tensor (including 
z turbulent diffusivity) . 

H(z) is the source of heat per unit volume. 

a, • dA 
A ,s Hi* 

I is the bottom- surface source of mass p«r unit area. 

Q is the bottom-surface source of heat per unit area. 

The following variables occur in the governing equations and are 
defined thus : 

Density, a, is assumed to vary with temperature in the following form: 

p = 1.02913 + (-2.0 x 10" 5 }T + (-4.8 x 10 -6 )T 2 gm /cc (4.5) 

Equation (4.5) and the discussion of the constants were given by Sengupta 
and Lick (1974). 

Eddy diffusivity, K , is a function of both thermal and current structure 
of a lake. The form used here was deduced by Rossby and Montgomery 

(1935). 

K z = K zo (1 + a 1 R i )l (4 * 61 

where, 

K is the eddy diffusivity under neutral conditions. 

Sundaram et al. (1969) gave an empirical form for the yearly variation 

of K . 
zo 



A sinusoidal variation of K with a mean 0.21 cm 1 /sec and amplitude of 
0.0S2 cnt 2 7sec is assumed in°Equation (4.7). 

The semi-empirical constant a, = 0.1 is estimated in this study by 
comparing the values used by Sundaram et al. (1969) and the original 
value used by Monin and Obukhov (1954). 

The third variable appearing in Equation (4.6) is R., the Richardson 
number. This is defined as: 

a v 9 * 2 3T 

where, 

a is the volumetric coefficient of expansion of water. 

a is defined as: 

a = 0.0 + 1.538 x 10" 5 (T-4) + (-2.037 x 10~ 7 )(T-4) 2 (4.9) 

The constants in the above equation were taken from Sundaram et al. 
(1970). a can also be estimated from Equation (4.5); however, the 
use of Equation (4.9) leads to a higher accuracy. 

g is the acceleration due to gravity. 

w* is the friction velocity calculated from known conditions above the 
lake. 

The heat source, H, is that part of the solar radiation transmitted 
exponentially through the depth of the lake. The relation used in this 
study has been used by Dake and Harleman (1969) and by Mitry and 
Ozisik (1976) : 



H(z) = n(l-S)AU) * EXP(-n(z-h)) (4.10) 

o 



where, 



8 - fraction of the solar radiation absorbed at the surface, (0.5). 

n = extinction coefficient, (0.75 cm ). 

h = depth of the lake, (m). 

i = net solar radiation reaching the water surface, (cal/cm 2 /day) . 

Equations (4.3) and (4.4) are parabolic and mathematically repre- 
sent a diffusion process. The solution of these equations requires one 
initial condition and two boundary conditions. 



Initial Conditions 

The temperature of the lake at spring homothermy is taken as the 
initial temperature. 

Boundary Conditions 

Surface: 

q = pC K |^, . = K (T e - T ) (H-.11) 

H s p p z 3z I z=h s E s 

where, 

T- = equilibrium temperature, (°C). 

T = surface temperature, (°C). 

K = surface heat exchange coefficient, (cal/cm 2 -s) 

T_ Evaluation 

Brady et al. (1969) showed empirically that fluctuations in the 
equilibrium temperature may be conveniently estimated using the approxi- 
mate relationship: 

T c = T . + H /K (4.12) 

c a s s 

where, 

T . is the dewpoint temperature, estimated from known conditions above 
the lake. 

H is the qross rate of short wave solar radiation. 

s a 

Since the dewpoint temperature tends to remain relatively constant through 
a single day. Equation (4.12) indicates that the main source of hourly 
fluctuations in T p is the solar radiation component. This generally 
reaches a maximum at solar noon, unless variable cloudiness interferes. 
At nighttime, T_ approaches the dewpoint temperature, which acts like 
a relatively invariant datum for periods of 24 hours or less. On an 
annual basis, *• <wever, both T . and H are generally much greater in 
summer than in winter. The dominant 'contribution to the amplitude of 
seasonal fluctuation in T p is the dewpoint temperature. 

K Evaluation 

The form used by Edinger and Ceyer (1967) was also used in this 
study: 

K s = 4. 5 + 0.05 T g + 3f(w) + 0.47f(w) (4.13) 



3 is found by applying standard curve-fitting techniques to pub- 
lished data pertaining to saturated vapor pressures at various tempera- 
tures; a convenient representation given by Edinger and Ceyer (1967) 
is: 

= 0.35 + 0.015 T + 0.0012 T 7 (mmHgy°C) (4.U) 

m m 

T <+ T H 

m 2 

The evaporative wind speed function f(w) used is also similar to that of 
Edinger and Ceyer (1967). 

f(w) = 9.2 + 0.46W 2 (W ~ 2 mmH "*) (4.16) 

m g 

where W is the wind speed (m/s) 

Bottom : 

The second boundary condition is at the bottom surface of the lake 
which is assumed to be perfectly insulated. 

7^- ! n = C»- 17) 

3r ' z=0 

Numerical Method 

Numerical integration of the governing equations require that these 
equations be replaced by finite difference equations. 

In this study, the following notations are used in the finite differ- 
encing of the basic equations: 

M The thickness of a slice (the basin is divided into equal number 

of slices). This is used as depth increment. 

th 

.it The time increment chosen at the n inteqration step. 

n ^ 

n, n-1 The present and last time at which computations were performed, 
j The space subscript. This is used as 

T. n = T(jAz, t n ) = T(z, t) 
where z is the vertical coordinate. 



Equation (4. 4) can now be written in difference form as: 
_ n _ — n-1 



where, 



- kb""' ih "» atv ,'m -'^.'h 1 



+ '^7' n "'23T"^v i ) ) n -'-t»AV z ,- , i 

Q, . , 






(pAK^T)! 1 ^ 1 = pAK r [(jH)Az, t^HTJ 1 ', 1 -Tj 1 " 1 ] 
(pAK^-nj^ 1 = pAK z [(j-})A2, t^.jHT?" 1 - TJ^ 1 ] 

9 z j - 



K = K [1 + aA- ^i)" 1 (a ).z. I^- 1 

Z. , , 20 1 W* 2 VII 3 z 

j+1 ' ' 

= K zo t(j^)iz]{1 + oA- ^^V l'^tV-tm.' 1 +T n + - 1 )](iHAz} 



(a v } j = A 1 +A 2 (T j " *> +A 3 (T j ■ 4)I 



<H = B i +B 2 T H +B 3< T H )2 



A * 2 

A. , = -T- 

j-i h 



= (i-i)AiI^] 



Terms like AT. , and AT. . can be defined as 
1+1 j-i 
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iT W ■ T M " T i 



iT H ' T i " T i-t 



The sequence in which calculations are performed is as follows: 

The dependent variables, T, K , w*, a , p, T-, and K are 
initialized. The area of each slice is calculated. The title or the pre- 
sent year is printed. The values of the variables, K , w*, a , p, T_, 
and K , are then calculated. The temperatures of th* slices are finally 
calculated. If the temperature profile is unstable, mixing of the unstable 
portion of the profile is undertaken. 

During the next time step, the temperatures are updated, and the 
dependent variables are calculated again. 

The values of the temperature, T, eddy diffusivity, K , number of 
days and the surface heat transfer coefficient, K , are printed every 
time step, every day or normally at the end of each month. At the end 
of the present year, the new-year title is printed and computations con- 
tinue as listed above. A detailed description of the programs and flow 
chart has been published in a user's manual, Lee and Sengupta (1980). 
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SECTION 5 
MODEL CALIBRATION 



The model described in the previous section was applied to Cayuga 
Lake, New York. The main reason was to calibrate the model. Cayuga 
Lake was chosen for this calibration because of the vast amount of nu- 
merical work done on this lake. These include Sundaram and Rehm (1971), 
Sundaram et al. (1971), Sundaram and Rehm (1972), and Mitry and Ozisik 
(1976). The model w*<s also verified using Lake Keowee data base; this 
is described in the nw <t chapter. 

CAYUCA LAKE APPLICATION 

Cayuga Lake is one of the finger lakes in west central New York. 
It is the longest, about 40 miles (64 km), of the group, and the second 
largest, about 66 square miles (171 sq. km), with steep banks cut by 
inflowing streams. It averages 2 miles (3 km) wide about 500 ft (150 
meters) deep, with a surface elevation of 381 ft (115 meters). Tang- 
hannock Creek and its 210-foot (64-meter) cataracts flow into the hke. 
The lake is drained to the north by the Seneca River, which joins the 
Oswego River and flows into Lake Ontario. Cayuga Lake is connected 
by canal with Seneca Lake on the west and Erie Canal on the north. 
The city of Ithaca lies at the southern end of the lake. 

INPUT QUANTITIES 

In the numerical integration of the governing equations fb>- Cayuga 
Lake, the following relations were used. The numerical values were 
taken from Sundaram et al. (1971) and Mitry and Ozisik (1976). 

Density, p, varied as in Equation (4.5), Section 4, of this report. 

Eddy diffusivity varied as in Equations (4.6) and (4.7) with 
a. =0.1 and the Richardson number as given in Equation (4.8). 

The volumetric coefficient of expansion of water, a , used is shown 
in Equation (4. 9). v 

A sinusoidal relation was assumed for the variation of the friction 
velocity, w*, with time; normally, w* is defined as 



w*=VfiJ7pT (5.1) 
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The form used is 

w* = A^ + B 1 »in(||| ♦ C } ) (5.2) 

where, 

A. = Average value of w*, (3.018 cm /sec). 

B- = Half the annual variation of w*, (0.762 cm /sec). 

C- = Phase angle, (2. 61) (chosen In such a way that at time t = 0, 



w* = initial value of the friction velocity). 



The heat source^.H, is taken to vary as in Equation (4.10), with 
S = 0. 5, n = 0. 75 cm and h = 70.0 m ( = 200 ft). 

The net solar radiation, 4 , reaching the water surface Is taken 
from Mitry and Ozisik (1976): ° 

* o =A 2 + B 2 sin(^t + C 2 ) (5.3) 



where, 



A 2 = Average value of 4 , (6. 14 x 10 cal/cm 2 s). 

-3 

B. = Half the annual variation of 4 , (3.52 x 10 cal/cm 2 s). 

C, = Phase angle, (0.049) (chosen in the same way as C.). 
The heat flux, Q, in Equation (4.4) is defined as 



Q = ( C C ATQ J/(A(z)) (5.4) 

P P 



where. 



Q = Volumetric discharge from the power plant. In this study 
p Q = 1.508 x 10 8 cm 1 /sec; this value is chosen to correspond 
to the pumpina velocity used by Sundaram et al. (1971), (1/4 
ft /day) . 

.iT = 10°C = Assumed temperature change through the condensers 
of t'^# power plant. 

The pumping velocity V is defined as 

- » 

The pumping velocity term effects are felt between the intake and the 
level at which the heated effluent becomes neutrally buoyant (effective 
discharge level). 
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_ a 

K * Assumed constant at 5. 65 x 10 cal/cm 2 -s. 

s 

-4 
For a postulated 3500 MW plant for Cayuga Lake, a 8. 79 x 10 

cal/cm 2 s of waste heat will have to be rejected. 

Surface area of Cayuga Lake is 66 sq. miles (171 sq. km). 

The intake to the power plant is fixed at 38.1 m (a 125 ft) from 
the surface of the lake. 

Initial temperature is taken as the temperature of the hke at spring 
homothermy (usually coincides with the minimum surface temperature 
attained by the lake). For Cayuga Lake, this occurs around March, 
and the temperature assumed is T = 2. 9°C. 

Finally, a sinusoidal variation is assumed for the variation of the 
equilibrium temperature with time. 

T E =A 3 *B 3 sin(J£t ^ (5.6) 

where, 

A. = Average value of T-, (11°C). 

B, = Half the annual variation of T_, (16°C). 

C = Phase angle, chosen in the same way as C 1 or C_. 

The justification for using Equation (5.6) is described in Appendix 
D of this report. 

Two topographies are studied in the application of this model to 
Cayuga Lake; these are a right circular cylinder and a circular para- 
boloid approximations for the lake. 

Cylindrical Topography 

The area (radius, B = 7. 38 x 10 cmLpf the lake is assumed con- 
stant throughout the depth. Term A' or -t— = 0; see Equations (4.3) 
and (4. U). 

Circular Paraboloid Topography 

The lake J$ assumed to be a circular paraboloid with surface radius, 
B = •*. 38 x 10 cm. The area at any depth, z 'measured from the deep- 
est punt of the take) is given by: 

A'=^- (5.7) 

Thus, A' is a constant: 
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A'=^ (5.8) 



RESULTS 



Computations for a yearly cycle for Cayuga Lake are presented. 
The verification data base consists of vertical temperature profiles com- 
piled oy Henson et al. (1961). The comparison of simulated and observed 
vertical temperature profiles are shown in Figures 1 , 2 and 3. Each 
figure shows five profiles representing observed (M), discharge (cylin- 
drical domain - CD); discharge (parabolic domain - PD) ; nc-discharge 
(cylindrical domain - CN) ; and no-discharge (parabolic domain - PN). 

The no-discharge simulations are in good agreement with the data 
(the data was for no-discharge conditions) . 

The parabolic case has somewhat better agreement since it represents, 
qualitatively, the decrease in area with depth. However, the closeness 
of the simulated results for the two cases is surprising. Most lakes have 
the rate of decrease in area with depth greater than a paraboloid, which 
has a linear decrease. Thus, when realistic area changes are used, a 
greater difference between cylindrical and paraboloid cases can be ex- 
pected. 

The effect of discharge is significant only in the top layers till July. 
This is because the heated discharge rises to the surface during the first 
half of the year. For the later months, the discharge temperature is 
lower than the surface temperature causing the discharge to reach static 
equilibrium somewhere below the surface. Thus, significant thermal 
effects of discharge are seen at mid-depths till December. The tempera- 
tures are higher at these depths for the paraboloid topography. !n 
general, a temperature difference of the order of 3°C over the no-dis- 
charge case, can be seen. At the end of the annual cycle, a residual 
temperature increase of 1.75°C is detectable. 

In the above paragraph, the discharge from the power plant is 
treated as a plane source which is injected into the domain at a level 
where the discharge temperature equals the temperature of that level. 
This level could be called the effective discharge level. The effects of 
the pumping velocity term, therefore, can only be felt from the intake 
level to the effective discharge level. A temperature rise of 10°C through 
the condensers was assumed This temperature difference justifies the 
use of density as a function of temperature. 

Figures t and 5 show the annua! stratification cycle. From these 
figures it is evident that the surface temperature difference between the 
four cases is less than 2°C over the annual cycle. However, at mio- 
depths the paraboloid discharge case shows a 5°C difference compared 
to the no-discharge case; for the cylindrical discharge case this difference 
is 3°C. Generally, the highest surface temperatures are reached after 
150 days, while the highest equilibrium temperature occurs after 120 days. 
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Thus, there is approximately a 30 day tag in the surface temperature 
response. The maximum temperatures at mid-depth occur after 210 days 
for the no-discharge case; for the discharge case, the corresponding time 
is 210 days. No significant phase lag between cylindrical and paraboloid 
cases are observed. 

Figures 6 and 7 show the eddy diffusivity variation with depth 
(time in days used as parameter} for the cylindrical and paraboloid 
cases respectively, 't is observed that thermal discharge causes in- 
crease in eddy diffusivity in the epilimnion owing to increased mixing. 
No significant changes are seen in the hypolimnion. There is a temporal 
increase in the differences between the discharge and no-discharge cases. 
These observations are true for both domains. However, the diffusivity 
values are larger for the paraboloid cases. Also, at any given time, 
the paraboloid case shows deeper thermoclines. 
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Figure 1. Vertical temperature profiles, Cayuqa Lake 
(from to 90 days) 
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Figure 2. Vertical temperature profiles, Cayuga Lake 
(from 120 to 210 days) 
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Figure 3. Vertical temperature profiles, Cayuga Lake 
(from 2U0 to 330 days) 
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Figure '4. Stratification cycle cylindrical domain. Cayuaa Lake 
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Figure 5. Stratification cycle parabolic domain, Cayuga Lake 
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Figure 6. Variation of eddy diffusivity with depth cylindrical domain, 
Cayuga Lake 
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Figure 7. Variation of eddy diffusivity with depth paraboloid domain, 
Cayuga Lake 
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SECTION 6 
MODEL VERIFICATION 



The model calibrated in the last section was applied to Lake Keowee. 
There were two objectives. The first was to simulate the stratification 
behavior of Lake Keowee for the period 1 977 to 1979, the entire existence 
of this river-dammed lake, which started receiving power plant heated 
discharges in 1973. The second objective was to test the accuracy of 
this model for a situation where horizontal uniformity (see Appendix A) 
of temperature was somewhat suspected since there are two basins con- 
nected by a canal. 

LAKE KEOWEE APPLICATION 

Lake Keowee is located 40 km west of Greenville, South Carolina. 
It is the source of cooling water for Oconee Nuclear Power Station (ONS). 
It was formed from 1968 through 1971 by damming Little and Keowee 
Rivers. A connecting canal (maximum depth 30.5 m) joins the two main 
arms of the lake. Flow out of the lake, except for negligible leakage 
through Little River dam, is through the Keowee Hydro Station. Lake 
Keowee also exchanges water with Lake J ocas see- pumped storage station. 
The three-unit steam electric ONS with a net capacity of 2580 Mwe 
started operating in July 1973. ONS operated at annual gross thermal 
capacity factors (CTCF) of 11, 28, 69 and 59% in the years 1973 through 
1976, respectively, Duke Power Company (1976). From 1977 through 
1979 the CTCF varied from 65 to 75 , Duke Power Company (1979). 

INPUT QUANTITIES 

Since long term effects are the important factors in this sutdy, the 
input data to the one-dimensional model were mainly monthly averages 
supplied by Duke Power Company (1976 and 1979). These are: 

1. Maximum depth of Lake Keowee = 45 m (si 50 ft). 

2. Initial temperature (minimum surface temperature of Lake Keowee) 
= 7.5°C. 

3. ONS condenser cooling water flowrate, Q , Table 1. 

4. Temperature rise in condenser cooling water, Table 2. 

5. Outflow through Keowee Hydro Station, Table 3. 
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6. Flow through Jocassee-pumpted storage station, Table 4. 

7. Wind speed, Table 5 (0.2 meters from surface). 

8. Solar radiation reaching the water surface, $ , Table 6. 

o 

9. Dewpoint temperature. Table 7. 

RESULTS 

Computations were carried out for nine years, 1971 through 1979, 
using the above monthly-averaged data. The measured temperature 
profiles referred to in this discussion are horizontal -averaged data col- 
lected from seven different locations in Lake Keowee, as shown in Figure 
8. This is in approximate agreement with the lateral uniformity assumed 
in the formulation of this and other one-dimensional models. 

Computed (with discharge) and measured temperature profiles for 
1971 to 1979 are presented together whenever possible, Figures 9 through 
17. A second set of figures, 18 to 26, for the nine simulated years show 
the computed t .mperature profiles (with and without discharge) . Other 
figures presented include: a typical variation of eddy dlffusivity with 
depth (time in days as parameter). Figure 27; and a typical temporal 
variation of the depth of the thermocline. Figure 28. Finally, the 
variation of the surface heat exchange coefficient with time in years 
(Figure 29), and the stratification cycle for Lake Keowee (1971-1979), 
Figure 30, are presented. All the above figures are computer drawn. 

For 1971, the year the lake came into existence. Figure 9 shows 
the computed temperature profiles (solid lines) and the measured tem- 
perature profiles (broken lines). The surface temperature comparisons 
between the measured and computed are excellent except for the start 
of the heating season when the computed temperatures seem to lag the 
measured. The maximum lag occurs in April and reduces with the 
heating season and completely disappears in October when cooling is in 
full swing. The hypolimnetic temperatures show good agreement with 
measured data. The mid-layer temperatures compare perfectly for the 
first half of the year. From August to November, the trend changes 
and the computed profiles lead the measured. However, the modeling 
of the formation, maintenance and decay of the thermocline, is good. 
This can be seen by comparing the shapes of the measured and computed 
temperature profiles. The comparison improves with increasing years, 
Figures 10 to 17, and generally follows the same trend as outlined above. 
The comparison of measured and computed temperature profiles cannot 
be expected to be more precise for the following reason. The measured 
values were taken on some specific day (different for some months) 
during the month while the computed results show the temperature pro- 
files at the end of the month. 

To estimate the impact of ONS operations on the thermal dynamics 
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of Lake Keowee, a sat of simulations (with and without discharge) are 
presented. These are summarized in Figures 18 to 26. Two years are 
selected for detailed discussion, 1975 and 1978. The earlier year, 
Figure 22, was a year of very high ONS annual gross thermal capacity 
factor and the later year, Figure 25, is representative of the situation 
towards the end of computations. In Figure 22, although the discharge 
surfaced only during the early and later months of the year, the effect 
of discharge (broken lines) is clearly evident during most of the months 
at the subsurface layers (especially below the thermcline) . Figure 25 
shows the same situation for 1978. This time, however, the discharge 
surfaced from January through May and then from October through 
December, at total of three months more than in 1975. The hypolimnetic 
temperatures were also higher. Finally, by the end of computations, 
1979, there was a residual temperature difference of almost 3°C due to 
discharge. ONS started operations in July 1973, and accordingly, the 
first three years. Figures 18 to 20, show no differences between dis- 
charge and no-discharge simulations. 

Figure 27 shows the variation of the eddy diffusivity with depth 
for 1975, one of the nine years simulated. In this figure, tne maximum 
diffusivity always occurs at the surface for every month. The non- 
stratified months show a constant diffusivity while the stratified months 
show very high mixing in the epilimnion and very low mixing in the 
hypolimnion. 

Figure 28 shows the temporal variation of the depth of the thermo- 
cline for 1975. This figure is in agreement with one of the generally 
accepted theories of the formation, maintenance and decay of the thermo- 
cline. 

Finally, a summary of the stratification cycle for Lake Keowee for 
the nine years simulated is presented in Figure 30, and the variation of 
the surface exchange coefficient for the same period is shown in Fiqures 
29 and 30. The maximum surface temperatures occur around 240 days, 
and the minimum surface temperatures around 60 days for each year. 
The equilibrium temperature of the lake is superimposed. The maximum 
equilibrium temperatures occur around 210 days. The surface temperature 
therefore lags the equilibrium temperature by 30 days. Two curves are 
shown for the mid-layer temperatures to highlight the effect of ONS 
operations. The temporal increase of this effect can be seen clearly. 
The maximum difference occurs during the last quarter of each year. 
The variation of the surface heat exchange coefficient, as expected, 
follows the same trend as the variation of the surface temperatures. 
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TABLE 1. Oconee Nuclear Station Condenser Cooling Water Fiowrata (m 3 /min) 



Month 


1973 


1974 


1975 


1976 


1977 


1978 


1979 


January 





3069.3 


461 2. 4 


6069.3 


S045. 8 


6176.7 


7207. 7 


February 





3069. 4 


3694. 9 


4440.2 


4985. 2 


6444. 6 


7319.9 


March 





2976. 9 


5456.8 


4874.3 


5113.5 


51 95. 7 


741 9. 5 


April 





2807. 3 


5570. 8 


4272.1 


6013.6 


4811.8 


7275. 8 


May 





2164.6 


6494. 3 


3970.7 


6302. 4 


4984. 2 


4189.1 


June 





4171.8 


6574. 2 


51 97. 6 


4385. 3 


5659. 9 


5381.2 


July 


1890.2 


5334.6 


7104.2 


5830.0 


5038.6 


7058. 8 


4753.3 


August 


1910.3 


4727. 1 


7510.1 


7248. 3 


5708. 9 


7914.9 





September 


2099. 5 


5961.4 


7201.6 


6785.4 


6964.0 


6557. 3 





October 


2232.5 


4953. 4 


6993.4 


5637. 8 


6754. 7 


7407.4 





November 


2170.7 


4202.1 


7467.1 


5809.2 


4697. 6 


6065. 1 





December 


3284.6 


5225.6 


6850.9 


491 4. 8 


5854. 6 


6503. 5 
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TABLE 2. Oconee Nuclear Station Condenser Tamperatura Rise. AT(°C) 



Month 


1973 


1974 


1975 


1976 


1977 


1978 


1979 


January 


— 


4.2 


6.3 


10.6 


12.5 


9.0 


10.3 


February 


— 


7.4 


4.8 


7.3 


11.4 


11.0 


10.4 


March 


— 


8.4 


6.2 


7.1 


10.4 


13.2 


9.6 


April 


— 


8.0 


8.3 


5.1 


11.4 


9.7 


9.9 


May 


— 


2.7 


8.8 


5.8 


9.4 


10.1 


8.2 


June 


— 


6.0 


8.8 


9.3 


8.4 


8.1 


7.1 


July 


5.3 


5.0 


8.3 


7.4 


7.4 


7.9 


5.0 


August 


4.6 


4.8 


7.8 


8.5 


5.0 


7.5 





September 


5.3 


5.8 


7.4 


8.0 


5.0 


7.6 


— 


October 


7. 3 


3.5 


7.7 


7.8 


3.8 


6.2 





November 


7.7 


7.9 


8.5 


6.7 


6.2 


8.4 





December 


4.1 


5.9 


9.4 


8.4 


7.9 


7.2 
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TABLE 3. Monthly Average Flowrates (m 3 /sec) - Lake Keowee Hydro Station 



Month 


1971 I 


1972 


1973 


1974 


1975 


1976 


1977 


1978 


January 


— 


17.5 


22.2 


71.2 


27.0 


39.5 


35.5 


30.1 


February 




32.9 


17.2 


118.1 


50.8 


16.2 


33.6 


22.6 


March 


— 


9.0 


33.6 


20.7 


82.7 


42.4 


25.4 


25.7 


April 


— 


3.0 


30.1 


16.7 


42.5 


40.7 


SO. 5 


18.7 


May 


— 


21.2 


31.6 


35.9 


63.8 


53.3 


i 

21.6 


52.3 


June 


— 


22.9 


55.7 


11.8 


16.1 


53.2 


27.5 


43.1 


July 


— 


14.3 


9.7 


13.1 


20.3 


33.2 


20.1 


12.2 


August 


— 


15.8 


16.1 


28.8 


39.9 


10.8 


55.2 


28.7 


September 


22.6 


20.2 


18.4 


17.0 


35.1 


16.0 


16.5 


32.1 


October 


14.9 


6.1 


3.3 


5.1 


49.8 


39.2 


7.5 


25.2 


November 


18.5 


30.8 


5.6 


23.0 


47.2 


27.2 


22.5 


13.4 


December 


8.8 


17,3 


19.3 


12. a 


44.1 


30.1 


20.5 


13.4 


Annual Average 


16.2* 


17.6 


21.9 


33.7 . 


43.3 . 


35.2 . 


27.2 


32.1 



* Average of September through December. 



29 



TABLE 4. Lake Jocassee Hydro Flows (cfs) 



Month 


1973 


1974 


1975 


1976 


1977 


1978 


1979 


January 


— 








505 


-608 


187 





February 


— 








509 


206 


54 





March 


— 








562 


1181 


516 





April 


— 








-79 


1046 


370 





May 


— 








555 


27 


1762.2 





June 


— 








1032 


489 


180.2 





July 


— 








112 


-831 








August 


— 





-163 


-747 


-460 


— 





September 


— 





554 


39 


101 








October 


— 





765 


304 


152 








November 


— 





615 


-251 


448 








December 


— 





301 


-915 


461 
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TABLE 5. Lake Keowee, Wind Speed (on /sec) 



Month 


1971 


1972 


1973 


1974 


1975 


1976 


1977 


1978 


1979 


January 


6.69 


6.69 


7.22 


5.8 


6.3 


7.4 


8.04 


7.9 


8.6 


February 


9.3 


9.26 


7.3 


5.8 


7.6 


8.5 


8.4 


6.8 


7.2 


March 


9.2 


9.20 


7.1 


7.7 


9.6 


7.9 


7.7 


7.6 


7.9 


April 


8.72 


8.7 


8.44 


8.7 


7.6 


7.6 


7.6 


7.6 


7.6 


May 


7. S 


7.53 


6.83 


6.8 


4.8 


7.3 


6.2 


6.7 


6.7 


June 


5.65 


7.95 


3.04 


6.96 


5.82 


6.4 


6.7 


4.7 


4.8 


July 


6. CIS 


6.64 


5.32 


5.2 


5.1 


5.9 


5.8 


5.7 





August 


5.75 


6.07 


5.1 


5.87 


5.4 


6.7 


5.4 


5.1 





September 


5.77 


5.47 


6.8 


6.74 


7.3 


7.13 


5.3 


5.7 





October 


7.02 


7.17 


7.1 


5.7 


7.7 


7.21 


7.2 


6.6 





November 


7.53 


7.13 


8.14 


7.2 


6.9 


7.27 


7.5 


5.8 





December 


8.3 


6.8 


5.6 


6.9 


7.2 


8.2 


7.2 


7.3 
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TABLE b. Lake Keowee Gross Solar Radiation (Langleys) 



Month 


1971 


1972 


1973 


1974 


1975 


1976 


1977 


1978 


1979 


January 


167.0 


176.0 


162.7 


191.4 


191.4 


209.8 


205.5 


227.0 


208.0 


February 


264.4 


257.6 


279.5 


226.9 


226.9 


310.9 


317.6 


308.0 


251.0 


March 


264.4 


352.5 


348.5 


326,'' 


1?6. 1 


338.6 


328,5 


408.0 


373.0 


April 


457.0 


448.0 


449.3 


397.7 


3*7. 7 


496.9 


427.3 


429.0 


479.0 


May 


480.5 


433.6 


449.5 


436.0 


436.0 


448.4 


473.0 


513.0 


513.0 


June 


478.0 


564.3 


507.7 


559.3 


559.3 


480.2 


543.3 


598.0 





July 


409.0 


493.8 


496.9 


459.5 


459.5 


288.3 


551.8 


568.0 





August 


428.2 


453.5 


391.6 


480.0 


480.0 


480.4 


423.9 


461.0 





September 


329,0 


386.3 


338.4 


339.2 


339.2 


345.1 


350.7 


385.0 





October 


261.3 


298.1 


341.7 


302.5 


302.5 


287.5 


286.6 


369.0 





November 


247.7 


220.9 


247.6 


231.1 


231.1 


237.5 


196.2 


232.0 





December 


147.7 


148.0 


154.0 


181.9 


181.9 


195.0 


178.2 


19'i.O 






32 



TABLE 7. Lake Keowee Dewpoint Temperature (°C) 



Month 


1971 


1972 


1973 


1974 


1975 


1976 


1977 


1978 


1979 


January 


3.0 


1.67 


1.0 


8.2 


3.0 


-1.0 


-6.6 


-2.8 


-3.33 


February 


0.9 


-2.22 


-1.0 


0.0 


3.5 


3.2 


-2.78 


-5.0 


0.0 


March 


6. 3 


1.11 


10.0 


6.3 


2.2 


3.9 


6.0 


1.2 


5.0 


April 


7.5 


6.6 


7.7 


10.7 


7.2 


11.2 


10.2 


9.6 


9.2 


May 


17.2 


11.11 


14.3 


17.2 


17.5 


14.0 


15.4 


14.0 


14.0 


June 


18.8 


13.13 


20.25 


17.8 


19.0 


18.3 


18.0 


19.4 


19.4 


July 


20.0 


18.77 


22.2 


21.0 


21.3 


19.8 


20.2 


20.8 


20.7 


August 


19.44 


22.22 


21.7 


21.0 


21.0 


18.0 


20.7 


20.8 


21.0 


September 


18.33 


18.8 


20.8 


17.5 


16.2 


15.4 


18.7 


15.5 


15.7 


October 


13.88 


11.5 


13.5 


10.2 


12.4 


8.2 


9.2 


9.3 


10.0 


November 


2.88 


5.9 


7.2 


6.0 


7.9 


1.0 


7.0 


9.0 


10.3 


December 


5.5 


4.0 


3.2 


3.8 


2.0 


-1.5 


0.4 


0.4 


0.9 
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Figure 8. Map of Lake Keowee 
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Fiqure 9. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1971 



a 

— -a 






uj^b.oo 



JAN. 



o 



20.00 

TEMP.ICJ 



o 
■ «=> 






FEB 



^^.OO 20.00 

a TEMP.(C) 



a 



x £ 



MARCH 



w'b'.oo 20-00 

a TEMP.fC) 



o 
■ a 



X S 



/ 



u/V.OO 



APRIL 



20.00 

TEMP.(C) 



o> 



o 

. — % 
i— 



X S 




MAY 



[ti^'.OO 20.00 



o 

— o 



U. 

X S 



u/^00 




TEMP- (C J 



Uj u.uu 20.00 

a TEMP.rC) 




fu'V.OD 20.00 



TEMP-ICI 



UJ 

a 



20.00 
TEMP.ICJ 



o 



u_ 







uj^'-OO 20.00 



IEMP.IC > 



o 



x £ 




OCT. 



uV^'.OO 20.00 



TEMP.I C J 



o 









00 
TEMP 



NOV 



20.00 
(CJ 



o 



X S 



ft 



DEC 



.00 
TEMP 



20.00 
IC) 



Figure 10. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1972 



o 



ll. 



r s 



a 

UJ 

f3 



^ 



00 

TEMP 



jnw 



20.00 

fCJ 



o 



o. 



Ii- 






^^.00 






FEB 



20.00 



TEMP.IC) 



o 






/ 






MRRCH £ 



Si^'.OO 20.00 



TEMP.fCJ 



ra 
»-° 

a 



APRIL 



.00 20.00 

TEMP.fCJ 



O 
— .o 

■ o 



I— 



xe> 




[u^-OO 20.00 

o TEMP.fC) 



o 







u/V-OO 20-00 

a TEMP.fC) 




i^.OO 20.00 

TEMP.fC) 



^.00 20.00 

TEMP.IC) 



o 

■ a 



tu 



^ 




u/V.0Q 20-00 

a TEMP.fC) 




uj^-OO 20.00 

a TEMP.fCJ 



o 
a. 






f 



NOV 



{tj^'-OO 20.00 



o 

■° 

a. 



X S 



DEC 



TEMP.fC) 



Sj^Voo 20 .00 



TEMP.fC) 



Figure 11. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1973 
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Fiqure 13. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1975 
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Figure 11. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1976 
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Figure 15. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1977 
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Figure 16. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1978 
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Figure 17. Predicted discharge temperature profiles. Lake Keowee. 1979 
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Figure 18. Predicted no-discharae temperature profiles. Lake Keowee, 1971 
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Figure 19. Predicted no discharge temperature profiles. Lake Keowee, 1972 
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Figure 20. Predicted discharge (borkert lines) and no-discharae temDerature orofiles 
Lake Keowee, 1973 
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Ficjure 21. Predicted discharge (broken lines) and no-discharqe temperature profiles. 
Lake Keowee, 1971 
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Figure 22. Predicted discharge (borken lines) and no-discharae temperature oro files. 
Lake Keowee, 1975 
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Fiqure 23. Predicted discharge (borken lines) and no-discharge temperature profiles. 
Lake Keowee, 1976 
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Figure 21. Predicted discharge (broken lines) and no-dlscharqe temperature profile? 
Lake Keowee, 1977 
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Figure 25. Predicted discharge (broken lines) and no-discharqe temperature profiles. 
Lake Keowee, 1978 
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Figure 26. Predicted discharge (broken lines) and no-discharge temperature profiles, 
take Keoweo, 1979 
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Figure 27. Variation of eddy diffusivity with depth, Lake Keowee, 1975 
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Lake Keowee, 1 975 
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Figure 30. Stratification cycle. Lake Keowee. 1971-1979 
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APPENDIX A 
DERIVATION OF THE MATHEMATICAL MODEL 

DERIVATION OF THE VARIABLE CROSS-SECTIONAL AREA MODEL 

Assuming lateral uniformity (horizontally uniform lake) , the isotherms 
coincide with the isopycnals. The equations for such a model are obtained 
from the basic balance equations of mass and heat: 

|£ = -7 • pv (A.I) 

j^(pC p T) = V • pC p K • 7 • pC p TV + H (A. 2) 

where, 

p = density. 

t = time. 

V = velocity of flow. 

C = heat capacity. 

T = temperature. 

K = heat diffusivity tensor (including turbulent diffusivity) . 

H = source of heat per unit volume. 

Most models have started with the above equations and neglected 
all the terms involving the two horizontal components (or dimensions) . 
The neglect of the horizontal divergence in such simple models is not 
always justified. There are at least two reasons for the existence of 
horizontal divergence in real lakes. 

1. The variation of horizontal cross-sectional area of the lake with depth. 

2. The existence of sources of heat and matter efflux on the bottom of 
the lake at depths above the deepest point. 
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The need to include these in the diffusion equations of lakes was 
already felt; Lerman and Stiller (1969), Outton and Bryson (1962) and 
Tzur (1973). Only Tzur (1973) formulated corrected diffusion equations. 

Another complication in applying a 1-D modaS to a lake is caused by 
seiches and internal waves. When these are present, the height of an 
isopycnal loses its meaning. According to Sweers (1968), averaging 
simultaneous temperature profiles at different points of a lake or taking 
profiles during a short interval of time in one or more points can lead 
to distortion of the shape of the thermocline. He recommends averaging 
the depths of isotherms instead. In a lake of constant depth this is 
probably the best available approximation but fails for a lake of varying 
depth. 

The effects of area change with depth are included by the following 
treatments of Equations (A.I) and (A. 2). Integrating Equation (A.l) 
over the volume of water below height h measured from the deepest point 
in the lake yields 

M dV = A • p VdV 

y3t y 

Using Gauss theorem on the right hand side yields 

/||dV = -p/fl • VdS (A. 3) 

v 3t s 

where S is a surface completely surrounding the volume V, hence 

dS = dC + dA 

where C is the surface area of the part of the bottom of the lake that 
is bounded by the contour at height z. As a subscript, z marks the 
vertical component of a vector. 

Using dV = Adz in Equation (A. 3) yields 

h. 

/i£Adz = -p/ft • VdC - p/* • VdA 

o 3t C A 



i.e. 



!^{Adi.-5,V n dC-.VA(h) (A 

O 



Integrating Equation (A. 2) over the volume of water below height h 
measured from the deepest point in the lake yeilds 
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/frpC TdV = /(7- P C R-VTJdV - /7-pC T?dV + / HdV 



Applying the divergence theorem to the first two terms on the right 
yields 

/A(z)-~pC Tdz = /fl. (pC R-7T)dS - /(VpC TtydS 

h 
+ / A(z)H(z)dz 
o ' 

Using dS = dA + dC yields 

h 



; A(z)~-pCTdz = / (pC R- v T) dA + / ( P C R- 7 T) dC 
o at p ft P Z «P " 



h h 

- / (pCT9)dA - / (pC T9) dC 
o p z o P n 



h 
+ / A(z)H(z)dz 
o 
i.e. 



h 3 h 



/A(z) — pC Tdz = pC A(h][K-7Tl + / p C [K- 7T] dC 
«p p z Q cc n 



-pC p A(h)TV z - /p c C c T c V n dC 

h 
+ / A(z)H(z)dz (A 5) 

o 

where, 

z = the vertical coordinate measured upward from the deepest point of 

the lake. As a subscript, it marks the vertical component of a vector. 

C - the surface area of the part of the bottom of the lake that is bounded 
by the contour at height z. 
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n = subscript that marks the component of a vector that is perpendicular 
to the lake-bottom positive upwards. 

A(z) = the horizontal cross-section of the lake at height z. 

As all the properties of the lake are constant over the horizontal 
isopycnals, the integration over the isopycnals becomes simple multipli- 
cation by their areas. 

Differentiating Equations (A. H) and (A. 5) with respect to the height, 
a set of 1-D equations are obtained from Equation (A. 4): 



at 3z z c n a z 



(A. 6) 



The last term in the above equation is the mass addition term {source 
term) from Equation (A. 5): 



a !t» c p t= 1t pC p aik -' t) *-» pC p atv * 



*U e C c (R.7T) n -p c C c T c V n l£ + AH 



In the above equation. 



p C (R • VT) = Conduction 
c c n 



Heat addition terms 



(A. 7) 



p C T V 
c c c n 



= Convection 



Because the horizontal gradients vanish, Equation (A. 7) can be simplified 
further by noting that 



(R- 7T) = K ~ 
z z3z 



(A. 8) 



Following a similar treatment used by Y. Tzur (1973), the slope of the 
bottom of lakes is usually small; it rarely exceeds a gradient of 0.1, 
Assuming thif- and considering an elemental area shown below, 

dA/l 



dZ 




dC/i 
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- a dz , 



tan 8 - l-n 5 gradient of the lake bottom 



Cos 9 = Costtan" 1 l^] (A. 9) 



also 



Cos e = 5T (A,10) 



From Equation (A. 9) and (A. 10) 



From Equation (A.1t) 



!£ = Cosltan" 1 tfy (A. 11) 

= Cos [arctan (gradient)] 

= n * average of Cos arctan 
gradient of the bottom surface 
at that depth. 



dA = ndC 



or 



i.e. 



dA dC 



dC = J. dA 

dz ~ n dz 

Since Cos arctan (0.1) = 0.995, n a - = 1. 

n 

Alternatively, n can be incorporated into the surface sources (equation 
(A. 7)). 

The 1-D equations can now be written in more concise and conveni- 
ent forms: 

A ft = -fz AcV z +,A ' (A ' ,2) 

and 

A ir<pC T) = fi(pC AK |J) - |-(pC ATV J + QA« + AH (A. 13) 

3t p 3 Z p Z3Z 3Z p Z 

where, 
dz 
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I = bottom-surface source of mass per unit area. 

Q = bottom-surface sou.ce of heat per unit area. 

Boundary Contitions 

Surface: 

Equilibrium Temperature Concept 

The bulk temperature of a vertically-mixed body of water under 
natural conditions, Tn, tends to increase or decrease with time, accord- 
ing to whether the sum, EH, of its heat inputs (net solar and atmospheric 
radiation, thermal discharge and heat outputs; back-radiation, evaporation, 
and condensation) is positive or negative: 

dT 

-3-^ = EH/pC p h (A.TO 

Following a similar procedure used by Edinger and Ceyer (1967), Equa- 
tion (A. 14) can be transformed (without linearizing the temperature 
dependence of the components of ZH) to yield 

-af = K s (T E- T n> / P C p h (A ' 15) 

where, 

K = surface exchange coefficient that depends on the water temperature 
and wind speed 

T_ = equilibrium temperature, defined as the hypothetical water surface 
temperature at which the net rate of surface heat exchange would 
be zero. 

Equation (A. 15) expresses this definition of T_ for the particular case 
of a vertically-mixed water body 2*n = when T = T E ). Edinger and 
Ceyer (1967) show that the two parameters are intrinsically coupled 
together via water temperature and meteorological conditions. For this 
reason iterative techniques are often used when applying Equation (A. 15) 
in a predictive role and because the rate of convergence is extremely 
rapid. The first step often yields satisfactory results. 

Brady et al. (1969) have shown empirically that fluctuations in the 
equilibrium temperature may be conveniently estimated using the approxi- 
mate relationship: 

T E =T d + H s /K s (A. 16) 
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where, 

T. - dew point temperature. 

H = gross rate of shortwave solar radiation. 

Since the dewpoint temperature tends to remain relatively constant 
through a single day, Equation (A. 16) indicates that the main source of 
hourly fluctuations in Tp is the solar radiation components. This gene- 
rally reaches a maximum at solar noon, unless variable cloudiness inter- 
feres. At nighttime T. approaches the dewpoint temperature which acts 
like a relatively invariant datum for periods of 24 hours or less. On 
an annual basis, however, both T. H are generally much greater in 
summer than in winter, "he dormnan? contribution to the amplitude of 
seasonal fluctuation in T_ is the dewpoint temperature. 

Exchange Coefficient Evaluation 

The same form used by Edinger et al. (1967) was used: 



K = 1.5 + 0.05 T + Sf(w) + 0.47 f(w) (A. 17) 

5 m 



where. 



T = surface temperature. 

S is found by applying standard curve- fitting techniques to published 
data pertaining to saturated vapor pressures at various temperatures ; 
a convenient representation given by Edinger and Ceyer (1967) is 



where. 



S = 0.33 + 0.015 T + 0.0012 T 2 (mn»Hg/C) (A. 18) 

m m 






The evaporative windspeed function f(w) used is also similar to that 
of Edinger and Ceyer (1967): 

f(w) =9.2 + 0.46 WifWm^mmHg" 1 ) (A. 20) 

where, 

W = wind speed (m/S). 

The first boundary condition is. at the surface and can now be stated as: 
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Bottom : 

The second boundary condition is at the bottom of the lake which is 
assumed to be perfectly insulated 

•p-i -a a (A. 22) 

and tha initial condition is the temperature of the lake at spring homo- 
thermy. 

T lnlttol" T . (A - 23 ' 



67 



APPENDIX B 

MEASURED TEMPERATURE PROFILES 
(Individual Stations) 



68 






a 
— o 



o. 



»- 



X? 



jrn 



(^J.OO 20. 00 

o TEMP.(C) 






2cT.0U 



UMP.<L) 



CD 

—•a 



I 



^2 
•-*? 

Vb.uo 



UJ 



O 

— . c» 

■*? 



itnp 



£0-00 
fCI 



MRY Z 



£i 



ac boo 



JUNE 






— p- 
*0 



ou 
ltMP-CLi 






FtP. Z 



o 

— C 



hRRCH £ 



a 






20.00 

ItMP.fU 



JULY 



*§ 



a ltMP-'CI 



x=» 

iZ3 



4» 



BPRll 



Jh^.OO 20-00 

a iEnp.ru 



a 

—a 
■ o 



12 



t 
i 



RUG- 



o UHP.lCl 



CD 



»- 



I-*? 

1>J u 



^' 



4 
I 
I 



U0 

TtuP 



iL-OU 



SEPT. Z 









00 

ttfiP 



?000 

ICi 



o 

— o 

■ a 



OCT. Z 



NOV 



S,^ 00 40. 00 

a UMP-fC) 



OtL 






Jh^-OU 20- OO 

o unP-tCi 



Fiqure B-l. Lake K«sowee measured temperature profiles. 1971 - Station 500 
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Figure B-2. Lake Keowee measured temperature profiles, 1972 - Station 500 
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Figure B-3. Lake Keowee measured temperature profiles, 1973 - Station 500 
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Figure B-U. Lake Keowee measured temperature profiles, 1974 - Station 500 
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Figure B-5. Lake Keowee measured temperature profiles, 1975 - Station 500 
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Figure B-6. Lake Keowee measured temperature profiles. 1977 - Station 500 
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Ftqure B-7. Lake Kaowee measured temperature profiles, 1977 - Station 500 
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Fiqure B-8. Lake Keowee measured temperature profiles, 1978 - Station 500 
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Figure B-9. Lake Keowee measured temperature profiles, 1971 - Station 501 
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Figure B-10. Lake Keowee measured temperature profiles, 1972 - Station 501 



o 


t 


— o 


1 


o. 




*—* 




►- 




u. 


• 


1-° 





o 



JRN 



(Ij^'.OO 20.00 

a 1EHP.CC) 



i 6 ? 






n_ 



FEB 



_ c tl.O0 2000 

a lEftP.IC) 



o 

— Q 

-a 






MRRCH 



E/^j.oo 20 oo 

o 1EMP.IC) 



o 

— o 



»_o 



APRIL 



Si^'.oa 20.00 



TEHP.(C) 



^4 
U> 



o 



E§ 



MRY 



£,^.00 20- 00 

a 1EMP.CC) 



a 



o. 



i_ca 



Uj^j'-OO 20.00 

a TEMP-(C) 



JUNE £ 



/ 
/ 
i 
i 



x 2 

■ o 

Si^'.OO 20.00 

a lEKP-(C) 



T" 



o 

— O 



JULY LZ 



X° 



y 



RUG. 



£} c b^oo zb.oo 

a TEMP.(C) 



o 

— o 



o. 



-j-o 



o TEMP 



SEPT 



20-00 
(C) 



■ Q 

a. 



xo 



Vb.OO 



^ 



UJ 

o 



OCT. 



2U.0U 
1EMP.CC ) 



a 
o. 



5§ 



^ y NOV. 



i^'.OO 20-00 

a TEMP. CO 



a 
— o 



a. 






DEC 



Srb-oo 



20. 00 
1EMP-CC) 



Figure B-11. Lake Keowee measured temperature profiles. 1973- Station 501 
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Figure B-12. Lake Keowee measured temperature profiles, 1974 - Station 501 
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Fiqure B-13. Lake Keowee measured temperature profiles, 1975 - Station 501 
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Flqure B-14. Lake Keowee measured temperature profiles, 1976 - Station 501 
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Figure B-15. Lake Keov measured temperature profiles, 1977 - Station 501 
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Figure B-16. Lake Keowee measured temperature profiles, 1978 - Station 501 
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Figure 17. Measured temperature profiles. Lake Keowee, Station 502, 1971 



—a 

m O 



x £ 



20-00 



[L^'.oo 

Q TEMP. (C J 

Q 



oo 
en 






Oo 






Mfly 



_ r \l . 00 20 . 00 

« TEMP.fC) 



o 

—a 

Q. 



JAN. £ 



iS 



FEB 



a_ ^ 

uj U-00 20.00 

a TEMP.fC) 



Q 

— Q 
O 



X £ 



Sj^'.OO 20.00 



TEMP.rC) 



o 

o. 



i- 



XO 



a 



O 
— o 



MARCH 



00 20.00 

TEMP.fC) 



JUNE £ 



X S 



JULY 



•00 20.00 

TEMP.fC) 



a 
— o 

o. 






APRIL 



u/^-OO 20.00 

a TEMP.(C) 



x 2 



r 
i 



( 

! AUG. 



uj'V-OO 20.00 

° TEMP.rCl 



— -O 

X S 



SEPT 



^•00 20.00 

ta TEMP.fC ) 



o 



X S 



OCT 



Sj'V.OO 20.00 



TEMP.fC) 



o. 



X S 



u/b-OO 20.00 

o TEMP.fC) 



o 



NOV. £ 



X S 



DEC 



O-CV. > 

uj U-00 20.00 

a TEMP.fC) 



Figure B-18. Lake Keowee measured temperature profiles, 1972 - Station 502 
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Figure B-19. Lake Keowee measured temperature profiles, 1973 - Station 502 
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Figure B-20. Lake Keowee measured temperature profiles, 1974 - Station 502 
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Figure B-21. Lake Keowee measured temperature profiles. 1975 - Station 502 
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Figure B-22. Lake Keowee measured temperature profiles, 1976 - Station 502 
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Figure B-23. Lake Keowee measured temperature profiles. 1977 - Station 502 



o 

■ Q 






o. 



JflN 



^.00 20-00 

o TEMP.fC) 



Q 
■ Q 



i— o 



FEB 



u/b.oo 20.00 

a TEflP.fC) 



o 

-— 

■ a 



u. 

x 2 

uj u 
a 



f 



a. 



MARCH Z 



X S 



r 



APRIL 



00 20.00 Si'V'OO 20.00 

TEMP.(C) o TEMP. (CI 



Kl 



1° 

a. 






MAY 



uj^JO 20.00 

a TEMP.fC) 



Q. 



y 



/ 



JUNE 






i*^ .ao 2uTno 

a TEMP.(C) 



a 
— o 



Q.Q 



/ 2 

JULY |Z 

12 



f 



/ 
j 

AUG. 



00 20. 00 Jb'V^OO 20.00 

TEMP.(C) a TEMP.(C) 



— o 
■ a 



X S 



fb'Voo 20.00 

a TEMP.(C) 



a 

— o 
m o 

a. 



SEPT. £ 






I 
OCT 



uj'V-OO 20.00 



TEMP. (CI 



o 

o. 



x 2 



a 



NOV. £ 

12 



DEC 



/V.OO 20-00 Sj^.OO 20.00 

TEMP.fC) a TEMP.(C) 



Figure B-2*. Lake Keowee measured temperature profiles, 1978 - Station 502 
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Figure B-25. Lake Keowee measured temperature profiles, 1971 - Station 503 
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Figure B-26. Lake Keowee measured temperature profiles, 1972 - Station 503 
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Figure B-27. Lake Keowee measured temperature profiles. 1973 - Station 503 
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Figure B-28. Lake Keowee measured temperature profiles, 1974 - Station 503 
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Figure B-29. Lake Keowee measured temperature profiles, 1975 - Station 503 
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Figure B-30. Lake Keowee measured temperature profiles, 1976 - Station 503 
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Figure B-31. Lake Keov/ee measured temperature profiles, 1977 - Station 503 
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Figure B-32. Lake Keowee measured temperature profiles, 1978 - Station 503 
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Figure B-33. Lake Koewee measured temperature profiles, 1971 - Station 504 



o 
— o 



lj_ 






ii'V.oo 20.00 



unp.rci 



o 

— Q 
Q. 



XO 



i^-OO 20- 00 

> TEMP.(C) 



a 



FEB. t 



X S 



MARCH 



.00 20.00 

TENP.CC) 



CD 



a. 



X S 

Q_ 



RPR1L 



^■00 20.00 

a TEMP.fC) 



KJ 



o 

— o 



o. 



X S 



Mflr 



i^.OO 20.00 

o TEMP.fC) 



o 



i*V. 00 20.00 

» TEMP.(C) 



o 

— o 



o. 



X S 



a 



o 



; JULY t 



see 



RUG. 



00 20.00 u'b.OO 20.00 

TEMP.fC) o TEMP.fC) 



O 




— o 




M Q 


r- 


Q. 


s 
/ 


~^ 


/ 


t— 


1 

» < 


u_ 


fc_^ 


1 


^E§ 





SEPT 



Cj^-OO 20.00 

a TEMP.fC) 



O 


1 


O 


/ 


O. 






/ 




1 




/ 




,' OCT. 


O 


i 


o 


• 



l.j^.OO 20.00 

^ TEMP.(CJ 



o. 



X S 



Si^.oo 



NOV 



20.00 

TEMP.fC) 



a 
— o 



DEC 



ui'V.OO 20-00 



TEMP. (CI 



Figure B-31 Lake Keowee measured temperature profiles, 1972 - Station 501 
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Figure B-35. Lake Keowee measured temperature profiles, 1973- Station 504 
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Figure B-36. Lake Keowee measured temperature profiles, 197* - Station 504 
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Figure B-37. Lake Keowee measured temperature profiles, 1975 - Station 504 
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Figure B-38. Lake Keowee measured temperature profiles. 1976 - Station 504 
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Figure B-39. Lake Keowee measured temperature profiles, 1977 - Station 504 
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Figure B-40. Lake Keowee measured temperature profiles, 1978 - Station 504 
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Figure B-11. Lake Keowee measured temperature profiles, 1971 - Station 505 



CD 
— O 



^1 



►- 






U/V-oo 20.00 

o TEMP.fC) 



—a 
■*? 



•JHN • j^ 



a 



FEB. 



00 70.00 

!EMP.(CJ 



■ o 
o. 



u. 






O 

■ *? 



HRRCH £ 
3-0 



.00 20.00 

TEMP. (CI 



/ 



I 



RPRIL 



Uj'TJ-OO 

a 



20.00 
TEMP.(C) 



o 

— o 

■° 

o 

u. 

-J-Q 



ujTl.OO 20.00 

a lErtP.fCJ 



o 

— o 



MRY £ 



5S 



uVoo 



?o.oo 
lEMP.(C) 



o 
■ 9 



JUNE £ 



a 



J 2. 



00 20 00 

TEMP.fC) 



a 
■ 9 



jult z 



xo 



( 



RUG. 



uj fl.OO 20.00 

a TEttP.(C) 



a 



o. 



t_0 






SEPT 



Jij'V.OO 20.00 

a TEMP.(C) 



— o 
■ a 



53 



UJ 

a 



u/^-OQ 20.00 



TEflP.fC) 



o 
— o 



OCT. £ 



X S 



00 20.00 

lEMP.fO 



o 
■ 9 



NOV. C 






OEC. 



.00 20.00 

unp.(C) 



Figure B-H2. Lake Keowee measured temperature profiles, 1972 - Station 505 
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Figure B-H3. Lake Keowee measured temperature profiles, 1973 - Station SOS 
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Figure B-44. Lake Keowee measured temperature profiles, 1974 - Station 505 
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Figure B-45. Lake Keowee measured temperature profiles, 1975 - Station 505 
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Figure B-16. Lake Keowee measured temperature profiles, 1976 - Station 505 
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Figure B-17. Lake Keowee measured temperature profiles, 1977 - Station 505 
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Figure B-4B. Lake Keowee measured temperature profiles, 1978 - Station 505 
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Figure B-49. Lake Keowee measured temperature profiles, 1971 - Station 506 
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Figure B-50. Lake Keowee measured temperature profiles, 1972 - Station 506 
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Fiqure B-51. Lake Koewee measured temperature profiles, 1973 - Station 506 
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Figure B-52. Lake Keowee measured temperature profiles, 1974 - Station 506 
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Figure B-53. Lake Keowee measured temperature profiles, 1975 - Station 506 
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Figure B-5H. Lake Koewee measured temperature profiles, 1976 - Station S06 
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Figure B-55. Lake Keowee measured temperature profiles, 1977 - Station 506 
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Fiqure B-56. Lake Keowee measured temperature profiles, 1978 - Station 506 



APPENDIX C 

MEASURED TEMPERATURc PROFILES 
(Av«rag« of Stations 500 Through 506) 
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Figure C-1. Lake Keowee averaged measured temperature profiles. Stations 501-506, 1971 
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Figure 2. Lake Keowee averaged measured temperature profiles. Stations 501-506, 1972 



1 ' my "'****Q!$pW 



00 



CD 

■ o 



u_0 



JRN. 



iD^'. 00 zT) . 00 

a TEMP. I C) 



Q. 






( 



MRY 



uj^.OO 20.00 

a TEMP.(C) 






zb-oo 



u/V.oo 

a TEMP.(C) 

o 

■° 



w'V.OQ, 20-00 

a TEMP.(C) 



a. 



FEB. £ 



3:9 



a 






o 



MRRCH \Z 



ro 



rpril 



IK) 20.00 {b'V-OO 2b. 00 

TEMP.(C) a TEMP.fC) 



JUNE tZ 






J ° 

a. 

JULY £ 
xs 



./ 



/" 



RUG. 



00 20.00 Sj^'OO 20.00 

TEMP.(C) a TEhP.(C) 



CD 
.O 

o. 



10 



Si^.OO 20.00 

a TEMP.fC) 



o 

t-iQ 

o. 



SEPT. £ 



xo 



' OCT • 



Sj* = tj 1 7oo zTT.oo 



TEMP.(Ci 



o 



a. 






a 



m 2b. 00 
TEMP.(C) 






( Nov. XL 



xg 



DEC. 



Si c b7oo 2b. 00 
a TEhP.ICJ 



Figure C- 3. Lake Keowee averaged measured temperature profiles. Stations 501-506, 1973 
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Fiqure C-1. Lake Keowee measured temperature profiles, Stations 501-506, 1974 
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Figure C-5. Lake Keowee measured temperature profiles. Stations 501-506, 1975 
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Figure C-6. Lake Keowee measured temperature profiles. Stations 501-506, 1976 
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Figure C-7. Lake Keowee measured temperature profiles. Stations 501-506, 1977 
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Figure C-8. Lake Keowee measured temperature profiles. Stations 501-506, 1978 



